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Abstract 

We propose a non-adiabatic scheme to detect the the geometric phase for the macroscopic 
state with the Josephson junction system. After extending the scheme to the two qubit cases, 
we provide a type of non-adiabatic geometric C-NOT gate which plays an important role in 
quantum computation. 

Testing quantum mechanics laws in macroscopic level is an important and interesting topic jj], B. Geometric 
phase H ^, |^] plays an important role in quantum interferometry and many other disciplines. It should be 
interesting to detect the geometric phase in macroscopic quantum state. Recently, after the proposal of geometric 
C-NOT gateQ with NMR, there is a proposal^ to adiabatically detect the Berry phase with Josephson junction 
of asymmetric SQUID, and also, to make a fault tolerate C-NOT gate through condition Berry phase shift]?]]. 
However, one should overcome two drawbacks in the suggestions in doing the geometric quantum computation 
or detecting the Berry phase only. One is the adiabatic condition which makes such gate not practical, the other 
is the extra operation to eliminate the dynamic phase. This extra operation seriously weakens the fault tolerate 
property. In this letter, we give a simple scheme to detect the geometric phase or realize the C-NOT gate non- 
adiabatically on the dynamic phase evolution curve with symmetric SQUID. By this scheme, the detection 
or the realization can be done faster and more easily. In our scheme, the above mentioned two drawbacks of 
previous suggestions are removed. We believe our scheme has led the idea of geometric C-NOT gate much more 
practical than before. 

For the two level system, Geometric phase is equal to half of the solid angle subtended by the area in the 
Bloch sphere enclosed by the closed evolution loop of the eigenstate. Recently ||, it was proposed to make a 
fault tolerate C-NOT gate using Berry's phase, i.e., adiabatic and cyclic geometric phase. It is well known|| 
that, with single qubit rotation operation, C-NOT gate can be used to carry out all computation tasks. The 
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experimental results with systematic errors were obtained on NMR||. Very recently, there is a new proposal 
to test the geometric phase and the conditional geometric phase shift with asymmetric SQUID. This is rather 
significant because testing the geometric phase on macroscopic quantum state itself is important in quantum 
mechanics, and also, a C-NOT gate with Josephson junction is undoubtedly a real quantum C-NOT gate. 
However, both proposals^, [?j rely on the adiabatic operations. This is a bit unrealistic because the macroscopic 
state of Josephson junction dephases fastjl|. The dephasing effects may distort the laboratory observation of 
the Berry phase seriously. On the other hand, for the purpose of quantum computation, any quantum gate 
should be run as fast as possible. Besides the adiabatic condition, both of the previous proposals take extra 
operation to eliminate the dynamic phase. This extra operation is unwanted for a fault tolerate C-NOT gate 
because if we can not eliminate the dynamic phase exactly, the fault tolerate property is weakened. For these 
reasons one is tempted to set up a new scheme which does not rely on the adiabatic condition and which does 
not involve any dynamic phase in the whole process. 

Geometric phase also exists in non-adiabatic process. It was shown by Aharonov and Anandan^ that 
the geometric phase is only dependent on the area enclosed by the loop of the state on the Bloch sphere. In 
non-adiabatic case, the path of the state evolution in general is different from the path of the parameters in the 
Hamiltonian. The external field need not always follow the evolution path of the state like that in adiabatic 
case. So it is possible to let the external field perpendicular to the evolution path instantaneously so that there 
is no dynamic phase involved in the whole process. 

Consider a superconducting electron box formed by a symmetric SQUID( see fig. |l|), pieced by a magnetic 
flux $ and with an applied gate voltage V x . The device is operated in the charging regime, i.e., the Josephson 
couplings Ejq are much smaller than the charging energy E c h- Also a temperature much lower than the 
Josephson coupling is assumed. The Hamiltonian for this system isM p| [To[ 

H = E ch {n-n x ) 2 - Ej(<S>) cos X , 

Ej(<&) = 2i?,/ocos (tt-§^J, E c h is the charging energy and n x can be tuned by the applied voltage V x through 
V x = 2en x /C x ( see figure^). The phase difference across the junction \ and the cooper pain number n canonical 
conjugate variables n] — i. $o = h/2e is the quantum of flux. So here n x and $ can be tuned externally. 
As it was pointed out earlier |7j ^ when n x is around 1/2, only two charge eigenstates |0 > and |1 > are 
important. The effective Hamiltonian in the computational two dimensional Hubert space is H = — • a, 
where the fictitious field B = [Ej,0, E c h(l — 2n x )] and a are Pauli matrices. In this picture, the state of 
cooper pair (| >) or 1 cooper pair(|l >) is expressed in | J,> or | |> respectively, where states | |> and | |> are 
eigenstates of Pauli matrix a z . Suppose initially $ = n/2 and n x = 0, the initial state is ji/'o >= |0 >= | !>. 
To avoid the dynamic phase, we select the geodesic curves in the time evolution. We can use the following 
scheme (see fig. ||). 

1. We suddenly change n^Note that n x can be tuned through V x .) from to n = |(1 -5), 5 > 0. 
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2. Wait for a time 



T= - (1) 

^{E ch 5f + (2E Ja f 

After this the state is rotated along geodesic curve CBA for an angle 7T. See fig. ||. 

3. We suddenly change n x to \ + 5, and wait for a time r again. The state is rotated along geodesic curve ADC 
for an angle ir again. See fig. ^. 

4. Do the meassurement to detect the interference effects. 

In the above operation, the time evolution curve will enclose an area in Bloch sphere. The eigenstate basis 
|± > , which are eigenstates of Pauli matrix a y will acquire a geometric phase of ±7 = ±2tt9, respectively. If 
the initial state is |^>(0) >= |0 >, at time 2t it is 

|^(2t) >=U(2t)\0>, 

where U(2t) is the time evolution operator from time to 2t. It has the property 

U{2t)\± >=e ±l7 |± > . 

Using the relation 

|0>=--j=(|+>-|->) 

we have 

|V(2t) >= -^=U(2t)(\+ > +|- >) = -sin 7 | T> -cos 7 | |> (2) 

In general, suppose the initial state is ^(O) >= a\ |> +P\ |> with constraint of \a\ 2 + \/3\ 2 = 1, we have 

|*(2t) >= (acos7 - /?sin 7 )| |> - (a sin 7 + (3 cos 7) | [> . (3) 

We can see that, the non-zero 7 phase here will cause detectable interference effects. Specifically, if we choose 
appropriate 5 value so that 7 = ^, we can observe the flip to initial state. This property can be potentially 
useful to realize the fault tolerate NOT gate in geometric quantum computation (|, |tJ. Equivalently, we can 
take rotations to the state on Bloch sphere around x— axis for the same angle 9 instead of changing n x in the 
scheme. We do so in the following for the conditional geometric phase shift. 

So far we have a simple scheme to non-adiabatically detect the geometric phase for the macroscopic quan- 
tum state of symmetric SQUID. Scince it is operated non-adiabatically, the total time needed here should be 
comparable to one cycle in the experiment of ref|l]]. That is to say, much shorter than the decoherence time. 
We believe the experiment of our scheme can be done with the same set up in rcf . We have chosen the zero 
dynamic phase path|jll|. We need not have to make additional evolution loop to cancel the dynamic phase and 
preserve the geometric phase, which is obviously too 
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The above scheme for non-adiabatic detection of geometric phase to single qubit can be easily developed to 
the two qubit system to make a C-NOT gate with symmetric SQUID. Suppose we have two capasitively coupled 
symmetric SQUIDs( see fig[|) 1 and 2 with same Ejq and same E c h- For simplicity we call them as qubit 1 
and qubit 2, respectively. We will use subscripts 1 and 2 to indicate the corresponding qubits. For qubit 1( 
control bit) we set $i = 7r/2 and n Xi i = during the whole process. For qubit 2( target bit), we set $2 = 3>o/2 
and n X2 = 1/2 initially. In quantum computation, the state for qubit 1 and 2 can be either |0 >= | j> or 
|1 >= I J.>. The weak interacting Hamiltonian is 

Hi = A(m - n x ,i)(n 2 - n X)2 ) 

the constant A is the charging energy that derives from the capacitive coupling. We assume A is much smaller 
than E c h. In our setting, the Hamiltonian reads 



Through our setting, to qubit 2 the interacting Hamiltonian is an effective conditional Hamiltonian dependent 



qubit 1 is |0 >= I |>. 

We first rotate qubit 2 around x— axis for angle —8 see fig.( ||). Note that E c h is much larger than A, so 
the state of qubit 1 is (almost) not affected by any operation in qubit 2 in the whole process. The interaction 
Hamiltonian will create an evolution path on the geodesic circle ABC( see fig. ||). After time r = 7r/A, we 
rotate qubit 2 around x— axis for another angle — (it — 28). Again wait for a time r. Then rotate qubit 2 around 
x— axis for angle it — 8 to let the Bloch sphere back to the original one. After the above operation, if qubit 1 is 
on state | j>, an evolution path of ABCDA on the Bloch sphere is produced; if qubit 1 is on state | |>, nothing 
happens to qubit 2. 

In the scheme we have used the rotation operation to qubit 2, around x— axis. This can be done by suddenly 
changing $1 to and n x \ — -0- — h \ to qubit 2. The sudden change of n x 1 is to cancel the effects from the 
interaction Hamiltonian if qubit 1 is on state | j>, thus an exact —8 rotation around x— axis is guaranteed. 
Everytime after the rotation task around x— axis is completed, we always set the parameters of $1 and n x ^\ 
back to 7r/2 and 1/2, respectively. Note we need not worry about the case of | |> state for qubit 1. Because 
the interaction is in the case, qubit 2 is just rotated around an axis which is slightly different from x— axis 
and is then rotated back around the same axis. There is still no net effect to qubit 2 in the operation. 

Thus for qubit 2, the final state is changed by the unitary transformation U(2t) in the following way 



Hi = -Aci z • 02z- 



(4) 



on the specific state of qubit 1. Explicitly, it is 



iAcr 2 if state of qubit 1 is |1 >= | T>, and it is if state of 
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' cos7(^) i sin j(9) 



where M 



' i the geometric phase acquired for initial state \X, + > (point A in 

y i sin 7 (0) cos j(8) 

fig-D and 7 (0) = -20. We see |7(0)| = tt/4 makes a C-NOT gate here( see fig^). This corresponds to d = tt/8. 

Thus, our sceme can be used to make a geometric C-NOT gate, which is fault tolerate to certain types 
of errors(|, ^2). Since the adiabatic condition is removed, the total operation time needed here should be 
comparable to that of a normal C-NOT gate. Therefore we believe our scheme has led the idea of geometric 
quantum computation much closer to the practical application. The various parameters of a normal C-NOT 
gate are listed in ref|| ||, ||). It has been estimated there that the operation time can be much shorter than 
the decoherence time. Instead of the inductively coupled system, here we have used the capasitively coupled 
system, however, this is not an essential mordihcation. In the real experiment, one has to have certain readout 



facility. The single eletron transistor can be a good candidate! 13 . 
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Figure 1: SQUID with symmetric Josephson junctions. It consists of a superconducting box formed by 
a symmetric SQUID, pierced by the magnetic voltage V x is the applied voltage, which determines the offset 
charge n x . The device operates in the charge regime, i.e. Ejq « E c h- 
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Figure 2: A scheme for geometric phase detection on symmetric SQUID. By suddenly changing the 
paramecters of n x and $, we can get the sudden fictitious field in x — z plane in directions perpendicular to 
geodesic plane CBA and ADC, respectively. The fictitious field will create the evolution path of CBADC. The 
angle between geodesic plane ABC and the equator is 9. The solid angle subtended by the area of ABCDA is 
4tt<9 



6 




Ck 



Figure 3: The capasitively coupled SQUIDs. 
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Figure 4: Non-adiabatic conditional geometric phase shift acquired through dynamic phase 
evolution path. These are pictures for the time evolution on Bloch sphere of qubit 2 only in the case that 
qubit 1 is | f>. If qubit 1 is | |>, there is no net change to qubit 2. Picture a shows that after the Bloch 
sphere is rotated around a;— axis for —9 angle, the interaction Hamiltonian will rotate the Bloch sphere around 
z'— axis. At the time it completes 7r rotation, i.e. t = i/A we rotate the bloch sphere around x— axis again 
for an angle of — (n — 29), then we get picture b. In picture b the sphere is rotated around z" — axis by the 
interaction Hamiltonian. Note that point B in picture b has changed its position now. The geodesic cure CBA 
is not drawn in picture b. After time r we rotate the qubit 2 around x— axis for an angle of 7r — 9. 
Picture c shows the whole evolution path on the Bloch sphere, point A evolves along closed curve ABCDA, a 
geometric phase 7 = —29 is acquired. Point C evolves along the loop CFAEC, a geometric phase —7 = 29 is 
acquired. 
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